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Abstract 

We have studied a model of self-attracting walk proposed by Sapozh- 
nikov using Monte Carlo method. The mean square displacement 
(ii^(t)) ~ t'^'^ and the mean number of visited sites {S{t)) ~ t'' are 
calculated for one-, two- and three-dimensional lattice. In one dimen- 
sion, the walk shows diffusive behaviour with v = k = 1/2. However, 
in two and three dimension, we observed a non-universal behaviour, 
i.e., the exponent v varies continuously with the strength of the at- 
tracting interaction. 
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There is great interest in random walks and interacting walks [1,2, 3]. The 
mean square displacement of a random walk, follows a power- law 

~ i^^. The ordinary random walk(RW) is diffusive, with v — 1/2, in 
all dimensions. For a walk with repulsion, such as self-avoiding walk(SAW), 
the exponent v is greater than 1/2. Random walk on a fractal is anomalous 
with V < 1/2 [2, 3]. Various models of interacting walks have been studied, 
such as true self-avoiding walk[4,5], generalized true SAW[6], the Domb- Joyce 
model[7] , and an interacting walk with a weight factor p for each new site that 
the random walker visits[8,9]. A comparative study of interacting random 
walk models was performed by Duxbury ct.al. [10,11]. 

Recently, Sapozhikov proposed a generalized walk in which the probabil- 
ity for the walker to jump to a given site is proportional to p = exp{—nu), 
where n = 1 for the sites visited by the walker at least once and n — 
for other sites[12]. If u < 0, the walker is attracted to its own trajectory. 
This walk is called a self-attracting walk(SATW). Monte Carlo studies have 
suggested that z/ < 1/2 for -u = — 1 and u = —2 on two dimension and 
1/4 < p < 1/3 on three dimension [12]. However, Aarao Reis obtained 
non-universal behaviour of the SATW in one to four dimension using ex- 
act enumeration method[13]. Prasad et.al. concluded that the SATW in one 
dimension is diffusive [14]. 

In the present comment we report results of a Monte Carlo simulation for 
the self-attracting walk in one, two and three dimension. We find that SATW 
in one dimension is diffusive. However, SATW in two and three dimension 
shows non-universal behaviour. Monte Carlo simulations were performed on 
one dimensional lattice (10^-stcps with 2000-configurational averages for each 
parameter u), square lattice (10*^-stcps with 2000-averages), and cubic lattice 
(10^-steps with 2000-averages). The lattice sizes used in this simulation were 
L = 10^ (ID), 1024 X 1024 (2D) and 200 x 200 x 200(3D). We always used the 
periodic boundary conditions. We calculated the mean square displacement 
(i?2(t)) and the mean number of visited sites {S{t)). 

Fig. 1 (a) shows the log-log plot of the mean square displacement against 
the time. We obtained the exponent by a least-square fit: u = 0.500(7) 
for = 0, i/ = 0.500(9) for u = -0.5, u = 0.499(8) for u = -1.0, and 
v — 0.498(9) for u — —2.0. All the lines are parallel in the large-time limit. 
These results means that SATW in one dimension is diffusive and support the 
conclusions of Prasad et.al. [14]. This diffusive behaviour is further supported 
by the results of the mean number of visited sites, {S{t)), in Fig. 1 (b). 
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The mean number of visited sites {S{t)) shows the power law behaviour as 
{S{t)) ~ with k = 0.499(7) for u = 0.0, u = 0.499(4) for u = -0.5, 
k = 0.499(8) for u = -1.0, and k = 0.498(5) for u = -2.0. At short 
times the slopes of the mean square displacement depend on the parameter 
u. However, for large times all the lines are parallel and give the same slope 
regardless of the parameter u. Aarao Reis concluded that the exponent 
1/ in one dimension decreases continuously when u decreases, by the exact 
enumeration calculation up to t = 30 [13]. But their conclusion is not correct 
because the time steps are too short to reach the asymptotic behaviour, and 
the log-log plot of (i?^(t)) has curvature at short times. 

The log-log plot of the mean-square displacement and the mean number 
of visited sites of a two dimensional SATW versus time are shown in Fig. 
2 (a) and (b), respectively. The appearance of the plateau region for large 
times is due to the finite size of the substrate. In this region the walk touches 
the boundary. The exponents u were obtained as z/ = 0.500(3) for u = 0, 
u = 0.472(4) for u = -0.5, u = 0.404(5) for u = -1.0, and u = 0.300(9) for 
u = —2.0. Our results are consistent with those of Aarao Reis[13], Sapozh- 
nikov[12], and Lee [15]. Values of the exponents u decrease continuously when 
u decreases. We can not observe a critical value Uc proposed Sapozhnikov so 
that z/ = 1/2 for < M < 'Uc[12]. The mean number of visited sites {S{t)) 
shows the logarithmic correction for random walk {u = case) [1,2]. The 
slopes of the log-log plots of {S{t)) decrease continuously when u decreases. 

Fig. 3 (a) and (b) show a log-log plot for the mean square displacement 
and the mean number of visited sites for a three dimensional SATW in a cubic 
lattice. The plateau region is also due to the finite size of the substrate. The 
exponents u were obtained as z/ = 0.500(2) for u = 0, u = 0.485(4) for 
u = -0.5, ly = 0.466(5) for u = -1.0, and u = 0.294(6) for u = -2.0. The 
exponents u also decrease continuously when u decreases. These results are 
also consistent with those of Aarao Reis[13]. But they are not consistent 
with the prediction of Sapozhnikov that 1/4 < z/ < 1/3. In Sapozhnikov's 
argument there is an ambiguity regarding the scaling between the bulk cluster 
and the boundary cluster visited by the walk. We obtained the exponents 
k as k ^ 0.996(5) for = 0, A; = 0.993(2) for u = -0.5, k = 0.991(5) for 
u — —1.0, and k — 0.904(6) for u — —2.0. For random walk. A; = 1 in three 
dimension [1,2]. We observed that the values of k decrease slowly when u 
decreases. These Monte Carlo results are also in good agreement with the 
observations of Aarao Reis[13]. The continuous decrease of the exponent 
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V and k in three dimension is further confirmed by the simulation at the 
M-value between u = —1.0 and u = —3.0. We obtained the exponents as 
ly = 0.436(4), k = 0.970(3) for u = -1.5, ly = 0.240(5), k = 0.806(3) for 
u = -2.5, and ly = 0.190(8), k = 0.776(5) for u = -3.0. These results 
support that there are no crossover behaviours at the range —3.0 < u < 0. 

We calculated the exponents u and k in one, two and three dimensions 
by Monte Carlo simulations. We have concluded that SATW in one dimen- 
sion is diffusive. This observation supports the recent calculation of Prasad 
et.al.[14]. However, in two and three dimensions the exponents v decrease 
continuously when u decreases. These non-universal behaviours are consis- 
tent with those of Aarao Reis[13]. Our observations are limited on regular 
lattices. One still has to explore the scaling behaviour on general lattices 
such as fractal. 

This work has been supported by Inha University and by the Basic Science 
Institute Program, Ministry of Education, Project No. BSRI-97-2430. I 
wish to thank professor Vladimir Privman for his careful reading of this 
manuscript and good comments. 
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Figure Captions 

Figure 1: (a) Log-log plot of the mean square displacement versus time; 
(b) log-log plot of the mean number of visited sites versus time for one 
dimensional SATW ior u — (topmost curve), —0.5, —1.0, —2.0 (bottom 
curve) . 

Figure 2: (a) Log-log plot of the mean square displacement versus time; 
(b) log-log plot of the mean number of visited sites versus time for two 
dimensional SATW ior u — (topmost curve), —0.5, —1.0, —2.0 (bottom 
curve) . 

Figure 3: (a) Log-log plot of the mean square displacement versus time; (b) 
log-log plot of the mean number of visited sites versus time for three 
dimensional SATW ior u — (topmost curve), —0.5, —1.0, —2.0 (bottom 
curve) . 
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